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Abstract We have determined in detail the electron spectrum in the decay of bound 
muons. These results are especially relevant for the upcoming ji — e conversion ex- 
periments. 

1 Introduction 

From the observation of neutrino oscillations, we now know that lepton flavors are not 
conserved. However, the mixing and small neutrino mass differences seen in oscil- 
lations have a negligible effect on charged-lepton flavor violating (CLFV) reactions. 
Thus, the CLFV reactions provide a discovery window for interactions beyond Stan- 
dard Model expectations 1 1 2|. 

Muons play a central role in searches for CLFV |[Tl|2l, because they can be pro- 
duced in large numbers and live relatively long. One reaction that can be probed with 
particularly high sensitivity is coherent muon-electron conversion in a muonic atom, 

^l- + {A,Z)^e- + {A,Z), (1) 

where (A,Z) represents a nucleus of atomic number Z and mass number A. It has the 
advantage of producing just a single particle, a mono-energetic electron. It does not 
have the problem of accidental background that plagues searches for the decay jU+ — > 
e+7, which can be mimicked by a positron from a normal muon decay and a photon 
coming from the radiative decay of a different muon, bremsstrahlung, or positron 
annihilation-in-flight. Various experiments have been performed over the years to 
search for the conversion li3J. The most stringent results come from the SINDRUM 
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n Collaboration [4], which reports an upper limit of 7 x 10^'^ for the branching ratio 
of the conversion process relative to muon capture in gold. Several new efforts are 
being planned. In the nearest future, the DeeMe Collaboration |5| has proposed to 
reach 10^'^ sensitivity. Larger scale searches, Mu2e at Fermilab |6| and COMET at 
J-PARC (7], aim for sensitivities below 10^'^. In the long run, intensity upgrades at 
Fermilab and the proposal PRISM/PRIME at J-PARC may allow them to reach 10" 
sensitivity. A quite remarkable improvement of about four orders of magnitude, with 
respect to the current limit, is therefore envisaged. 

The success of the conversion searches depends critically on control of the back- 
ground events. The signal for the jj. — e conversion process in Eq. ([T]) is a mono- 
energetic electron with energy E^e, given by 

£^^1? — 'W^ — — /irec I (2) 

where is the muon mass, E\, ~ Z^a?m^/2 is the binding energy of the muonic 
atom, and E^ec — 'm^/(2?ma^) is the nuclear-recoil energy, with a the fine- structure 
constant and the nucleus mass. The main physics background for this signal 
comes from the so-called muon decay in orbit (DIO), a process in which the muon 
decays in the normal way, pL^ e^VeV^, while in the orbit of the atom. Whereas in 
a free muon decay, in order to conserve energy and three-momentum, the maximum 
electron energy is 111^/2, for DIO, the nucleus recoil can balance the electron's three- 
momentum taking basically no energy. This allows for the maximum electron energy 
to be E^e, close to the full muon mass m^. Therefore, the high-energy electrons from 
the muon decay in orbit constitute a background for conversion searches. 

2 Muon decay in orbit 

Several theoretical studies of the muon decay in orbit have been published. Expres- 
sions describing the electron spectrum including relativistic effects in the muon wave- 
function, the Coulomb interaction between the electron and the nucleus and a finite 
nuclear size have been available for some time ll8ll9l [T0l . However, the high-energy 
endpoint of the spectrum, which is the most important region for conversion-search 
experiments, was not studied in detail. References 111 1III2I did study the high-energy 
end of the electron spectrum, and presented approximate results which allow for a 
quick rough estimate of the muon decay in orbit contribution to the background in 
conversion experiments. We have performed a new evaluation of the DIO spectrum, 
considering in detail all the effects needed in the high-energy region ||T31 . Our re- 
sults describe the background contribution for — e conversion searches, as well as 
a check on previous low- and high-energy partial calculations 1 10 1 I j and an inter- 
polation between them. It is worth emphasizing that not only the high-energy region 
is relevant for conversion experiments, but the full spectrum is necessary in order to 
study reconstruction errors in the detector 

To obtain the correct result for the high-energy tail of the spectrum it is crucial 
to include nuclear-recoil effects, since they modify the endpoint energy (see Eq. (|2])). 
Also, to produce an on-shell electron with energy around m^, either the muon must 
be at the tail of the bound-state wavefunction or the produced electron must interact 
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Table 1 Values for muon energies E^, nucleai' masses nifj, and endpoint energies E^^. 



Nucleus 


Z 


£^ (MeV) 


iiiN (MeV ) 


E^,, (MeV) 


C 


6 


105.557 


11188 


105.06 


Al 


13 


105.194 


25133 


104.973 


Si 


14 


105.121 


26162 


104.91 


Ti 


22 


104.394 


44588 


104.272 



with the nucleus. This tells us that the full Dirac equation for the muon as well as 
the interaction of the outgoing electron with the field of the nucleus must be taken 
into account. Also, finite-nuclear-size effects will be most important in this region. 
Order a radiative corrections are not expected to significantly modify the results at 
the endpoint, and are not included in our results. Uncertainties in the modelling of 
finite nuclear-size effects induce errors in the spectrum that increase as we approach 
the endpoint, but those errors are never larger than a few percent. 

3 Results and discussion 

Here we present our results for the elements that are relevant for the upcoming 
conversion experiments. The DeeMe Collaboration plans to use a silicon-carbide 
target, whereas the Mu2e and COMET Collaborations are considering aluminum 
and titanium as targets. In Table [T] we give the values of the bound muon energy 
Efi = — £b and the electron endpoint energy Zs^g, for carbon, aluminum, sili- 
con and titanium. We present the results of the numerical evaluation of the spectra 
for those elements in Figs. [T] and |2] Regarding the nuclear distributions, a two- 
parameter Fermi distribution has been used for aluminum and titanium, a three- 
parameter Fermi distribution for silicon and a Fourier-Bessel expansion for carbon 
II13II14|[T51 . Those results are useful for assessing the DIO background events from 
carbon for the DeeMe experiment which will search primarily for conversion in sil- 
icon. They also illustrate the electron resolution requirements, as a function of stop- 
ping target, needed to reach future high sensitivity goals. 

The high sensitivity that the upcoming conversion experiments will reach may 
also allow them to improve the present bounds on some exotic muon decays, like the 
decay of the muon into an electron and a majoron (a Goldstone boson that appears in 
models where lepton number is a spontaneously broken global symmetry) l,16il . 
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Fig. 1 Electron spectrum, normalized to the free-muon decay rate Jo - The solid blue line is for carbon, the 
black dotted line for aluminum, the green dot-dashed line for silicon and the red dashed line for titanium. 
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Fig. 2 Electron spectrum for the full range of Eg ; see caption of Fig.^ 
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